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ABSTRACT 

In an early-type, massive star binary system, X-ray bright shocks result from the powerful collision 
of stellar winds driven by radiation pressure on spectral line transitions. We examine the influence 
of the X-rays from the wind-wind collision shocks on the radiative driving of the stellar winds using 
steady state models that include a parameterized line force with X-ray ionization dependence. Our 
primary result is that X-ray radiation from the shocks inhibits wind acceleration and can lead to a 
lower pre-shock velocity, and a correspondingly lower shocked plasma temperature, yet the intrinsic 
X-ray luminosity of the shocks, Lx remains largely unaltered, with the exception of a modest increase 
at small binary separations. Due to the feedback loop between the ionizing X-rays from the shocks 
and the wind-driving, we term this scenario as self regulated shocks. This effect is found to greatly 
increase the range of binary separations at which a wind-photosphere collision is likely to occur in 
systems where the momenta of the two winds are significantly different. Furthermore, the excessive 
levels of X-ray ionization close to the shocks completely suppresses the line force, and we suggest that 
this may render radiative braking less effective. Comparisons of model results against observations 
reveals reasonable agreement in terms of log(Lx/^boi) ■ The inclusion of self regulated shocks improves 
the match for kT values in roughly equal wind momenta systems, but there is a systematic offset for 
systems with unequal wind momenta (if considered to be a wind-photosphere collision). 

Subject headings: hydrodynamics - stars: winds, outflows, stars: early-type - stars: massive - X- 
rays: stars 
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1. INTRODUCTION 

A large fraction of massive stars reside in binary sys- 
tems, with r ecent estima t es of binarity for O-type stars 
of > 70% (IChini et all |20T1 ISana et al.l 12011) . In 
such systems, which consist of two hot luminous mas- 
sive stars, the collision of the powerful stellar winds 
leads to the formation of hi gh Mach number sh ocks that 
emit at X-ray wavelengths ([Stevens et al.lfl9"92T ). Histor- 
ically, colliding winds binary (CWB) systems have been 
characterized by high plasma temperatures and an X- 
ray over-luminos ity (compared to their expected single 
star brightness) (|Pollockl 119871 : IChlebowski k Garmanvl 
119911 ) with observ ational inferences corroborated by the- 
oretical models (iLuo et al.1 119901 : iStevens et al.l Il992t 
IPittard k Stevens! Il997fi ~ However, more recent stud- 
ies examining a wider population and using the XMM- 
Newton and Chandra satellites indicate that short pe- 
riod WR+O and O+O-star binary systems have a ra- 
tio of log(Lx/ £boi) — ~7, similar to that expected for 
single O-stars (lOwocki k Cohen] 119991: lDe~Becker et aT 
2001 IQskinoval I2005t ISana et al.l 120061: lAntokhin et al 
2008t iNazel 120091: iNaze et al.l 120111: IGagne et al.l l20ilT 



Gagne et al J 120121) . Therefore, superlative X-ray bright 
log(Lx/iboi) as high as -5 - appears to 
for the more massive CWBs with high 
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spread of over three orders of magnitud e in X-ray lumi- 
nosity from CWBs (|Gagne et al.l 12012] )? Specific stud- 
ies of archetypal systems around the higher luminos- 
ity end of the distribution have yielded promising re- 
sults. For example, three dimensional simulations of 
77 Carinae and WR140, which include orbital motion, ra- 
diative cooling, and in some cases radiative driving, are 
able to ex plain the X-ray lightcurves and spectra reason- 
ably well (lOkazaki et aTll2008t IParkin et al.|[200l 120111: 
Russell et alJ 120 111 ). In contrast, models of WR22 by 



Parkin k Gossetl (|2011fl over-predict Lx by up to two 



orders of magnitude, with the best agreement (a fac- 
tor of roughly six over-estimate) achieved when a wind- 
photosphere collision occurs and the majority of the X- 
ray emission is extinguished. Problems also arise for 
lower mas s CWB s ystem s. A model of an O6V+O6V 
binary by IPittard! (2 0091 ) and IPittard k Parkin! (j2010f ) 
revealed an estimated log(Lx/^boi) of between -6 and 
-6.3 (depending on the viewing angle). This should be 
compared against observed values for systems with or- 
bital periods of 2-3 days that have log(Lx/-^hni) between 
-6.2 a n d -7.3 (INazell2009l; IGagne et aill20lil IGagne et all 
120121) . IPittard k Parkin! (j201 01) have presented evidence 
that this discrepancy may, in part, be due to the spec- 
tral fitting procedure used to extract parameters from 
observations, which they show to under-predict the ac- 
tual X-ray luminosity (which is known from the models) 
by up to a factor of two, particularly for short period 
systems where occultation may occur. Alternatively, the 
inconsistency with observations may indicate that some 
additional physics is required in the models. 
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Fig. 1. — Cartoon illustration depicting the wind- wind collision in an equal winds massive star binary system without (left) and with 
self regulating shocks (right). The stars are represented by the circles. Arrows indicate the wind direction, and arrow length indicates 
velocity magnitude. The solid and dashed lines demark the regions of post-shock stellar wind and the contact discontinuity, respectively. 
The shaded region indicates plasma temperature - fainter shading corresponds to lower temperature. 



Consideration of radiative wind driving in a massive 
star binary system led to the discovery of two interest- 
ing effects: radiative inhibition ( Stevens fc Polloc9l994f) 
and sudden radiat ive braking (|Owocki fe Gavlevl 119951 : 
iGaylev et al.lll9 97). In the former, the acceleration of the 
stellar wind may be reduced by the radiation field of the 
binary companion, whereas the latter effect concerns hy- 
personic flows being effectively halted in their tracks en- 
abling a wind-wind collision in systems where one would 
not occur on the basis of a ram pressure balance alone. 
One factor that has not been previously studied in the 
CWB paradigm is the influence of the ionizing X-rays 
from the wind-wind collisi on shocks on the wind driving. 
iStevens fc Kallmanl (|1990l) examined the dependence of 
the radiative line force on X-ray ionization for the case of 
a high-mass X-ray binary system, finding that the stel- 
lar wind acceleration could be significantly suppressed by 
a particularly bright compact object because the exces- 
sive X-ray ionizati on reduces the radiative line force (see 
also !Stevenslll99ll ). This effect has also been explored 
for line-driv en instability shocks embedded in a massive 
star's wind (jKrticka fc Kubatl 120091: iKrticka et all 120091) 
and fo r radiatively drive n disk winds of active galactic 
nuclei (jProga et al.ll2000f ). 

In this paper we make the first attempt to examine 
the feedback of ionizing X-rays from the wind-wind colli- 
sion shocks on wind acceleration in a massive star binary 
system. Because of the direct coupling between the ra- 
diation force that drives the stellar winds and the ioniz- 
ing X-ray emission that results from the wind-wind col- 
lision, we term this effect self regulated shocks (SRSs). 
Fig. Q] depicts the basic scenario under consideration 
and highlights some key effects due to SRSs. Firstly, 
wind velocities are reduced (shorter arrows in the right 
panel) which causes a lower post-shock plasma temper- 
ature (fainter shading). Consequently, radiative cooling 
may become sufficiently important to intr oduce instabil- 
ities w hich will perturb t he sh o ck fronts (IStevens et alj 
fl99l IParkin fc Pittarl 12011 Ivan Marie et all 120111: 



iParkin et all I20TTI: lLamberts et~aTI |20TT[) . The goal of 
this work is to provide a qualitative picture, and initial 
quantitative estimates, of when/if the SRS effect might 
be important. Therefore, we will make simplifications in 
order to elucidate the physics. 

The structure of this paper is as follows: In § [2] we 
calculate the influence of X-ray ionization on the line 
force due to an ensemble of spectral lines. The semi- 
analytical wind acceleration model is described in § [31 
followed by results for model binary systems in § [4] An 
approximate model for SRSs is presented in § O We 
compare results to observations in § [5] and then discuss 
some implications of our findings, and possible avenues 
for going beyond the illustrative wind acceleration model 
adopted in this work, in §[7J The main conclusions of this 
work arc summarised in § [8] 

2. THE LINE FORCE 

For the wind models that will be presented in § [3j 
we need to compute the radiation force due to spec- 
tral lines for appropriate stellar parameters while ac- 
counting for the influence of X-ray irradiation arising 
from a wind collision. We will adopt an approximate 
treatm ent of the radiation force due to spectral li nes fol - 
lowing iCastor. Abbott, fc Klein (1975)1 (hereafter ICAKl) 
and c losely follow the approach by IStevens fc Kallmanl 
(|1990f ) to estimate the effect of X-ray ionization on the 
line force - essentially, our goal is to repeat their calcula- 
tions for the stellar parameters appropriate to our study. 
In this section we outline the method and the implemen- 
tation used here. For full details of the methodolo gy and 
discussions of its v alidity, we refer the reade r to ICAKl . 
lAbbottl (fl98l) and IStevens fc Kallmanl (fl99l . 

The total force due to lines is given by, 
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where o~ is the electron scattering opacity and F is the 
radiative flux. M(t) is known as the line force multi- 
plier, which depends on the dimensionlcss optical depth 
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parameter in a stellar wind, defined by 



' dv 

t = (7 C /Wth [ ~j- 



(2) 



where p is the mass density, w t h is the thermal velocity of 
a hydrogen atom and dv/dr is the radial velocity gradi- 
ent. The Sobolev optical depth of a spectral line between 
lower state I and upper state u is given by = r\ u jt 
where 



he niBi tU - n u B U} i 
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(3) 



Here, n\ and n u are the lower and upper level population 
number densities and B^ u and B u j are the usual Ein- 
stein coefficients for absorption and stimulated emission, 
respectively. The force multiplier, M(t), is composed 
from a sum over all line transitions 



M(t) = J2 a "d 



lines 



F v \- exp(-rj u ^t) 
F t 



(4) 



where Avd is the Dopplcr width and F v is the specific 
flux at the line frequency {y). 

To evaluate M(t), we need to supply a list of line tran- 
sitions (frequencies and oscillator strengths), specify the 
form of the radiation field F u , and compute the associ- 
ated level populations (n.;, n u , relative to the total den- 
sity p). _ 

The line list used in this study is drawn from two 
sources. For low-ionizati on metal atoms/ions, we use 
the CD23 line database of iKurucz fc Belli (j 19951) . From 
this source we include elements with atomic number 
6 < Z < 30 and include ionization stages I - v with the 
following exceptions: for C, we include only I - iv while 
for Z > 20 we include ions I - vn, where available. In 
order to extend our calculations to regimes of higher ion- 
ization, w e also included data fro m the CHIANTI atomic 
database (|Dere et al.lH997ll2009t) . From this source, we 
take line lists for H and He and the high ions of the as- 
trophysically abundance metals: C, N, O, Ne, Mg, Si, S, 
Ar, Ca, Fe and Ni (for each of these metal, we include 
CHIANTI line lists for all available ions that we did not 
take from Kurucz & Bell 1995; we excluded theoretically 
predicted lines from the database). In total, our line list 
contains ~ 7.7 x 10 5 transitions. 

The stellar radiation fie ld, F„ was taken from AT LAS9 
model atmosphere grids (|Castelli fc Kuruczl [2004') . For 
the specific stellar parameters used, see below. 

The level populations (n;, n u ) for each transition were 
comput ed in a two stage pr ocess. First we used Cloudy 
vlO.00 (jFerland et al.1 11998) to compute the ionization 
stage of a shell of gas illuminated by a specified radia- 
tion field. In all cases, we assumed that the irradiating 
spectrum contains two components: emission from the 
star and hard radiation associated with emission from the 
wind collision region. The shape of the stellar component 
was taken from the same model atmospheres used for 
F v j F . In setting the intensity of this component, we fol- 
low lStevens fc Kallman (1990) and consider only a single 
value for the ratio of the electron number density to the 
geometrical dilution factor (n e /W = 3.5 x 10 10 cm" 3 ). 



log I = -2 
log I = 
log I = 2 
log % = 4 




-10 



Fig. 2. — Force multiplier versus dimensionless optical depth pa- 
rameter (t) for log ^ = -2, 0, 2 and 4 (solid lines). Dashed lines 
show our fits (see text). These calculations are for a star with 
T eS = 38500 K, logg = 3.92 and log(Z/Z Q ) = 0. For compari- 
son, we also show calculations from Abbott (1982) for a star with 
T eff = 40000 K and logg = 4.0 (black spots). 



To describe the spectral shape of the hard ionizing radi- 
ation, we adopt a thermal Brcmsstrahlung spectrum at a 
temperature of 10 keV. The intensity of this component 
is specified as an ionization parameter, 



AirFxp-rnft 



(5) 



where, in this work, Fx is the flux of X-rays from the 
wind collision shocks, and p is the gas density. The value 
of £ is varied to quantify the affect of X-ray ionization on 
M(t). From the ion populations provided by the Cloudy 
calculations, we compute level populations assuming lo- 
cal thermodynamic equilibrium (LTE, adopting the gas 
temperature calculated by Cloudy). Although simplistic, 
this assumption makes it easy to compute the force mul- 
tiplier reasonably quickly. Ideally, full non-LTE calcula- 
tions should be performed for complete atomic models 
associated with each ion. This, however, would signifi- 
cantly complicate the calculation an d is not expected to 
quali tatively affect our findings (see [Stevens fc Kallmanl 
1990], for further discussion). 

2.1. Example calculation 

Using the procedure outlined above, we can calculate 
M(t) accounting for the effects of excess ionization (as 
controlled by £). As an example, we show results for a 
star with effective temperature T e g = 38500 K, surface 
gravity logg = 3.92 and solar metallicity log(Z/Z Q ) = 
in Fig. H 

As expected, our result s are in good agreement with 
iStevens fc Kallmanl (|1990f ). For each value of £, M is 
largest (and constant) when t is sufficiently small that all 
lines are optically thin. At large t, M decreases as lines 
become optically thick and saturate. For calculations 
with low ionization parameter (log£ < 0), M remains 
significant (M > 1) up to around t ~ 1 (M is essentially 
independ ent of ionization paramete r for log£ < —2). As 
found by IStevens fc Kallmanl ([l990h . we also see that as 
log £ is increased beyond zero, M drops and the regime in 
which M(t) is well-described by the optically thin limit 
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extends to higher t- values. For log£ > 3, the force mul- 
tiplier is always small. 

For comp arison, we also show in Fig. [2] the M- values 
reported by lAbbottl (|1982j ) from calculations for a star 
with similar parameters (T e g = 40000 K, \ogg = 4.00 
and n e /W = 1.8 x 10 11 cm~ 3 ). Since no ex cess ioniza- 
tion radiation was included by [Abbott] ([19821 ) , his calcu- 
lations should be compared to our results for the lowest 
ionization parameter shown (log£ = —2). In general, the 
agreement is very good - the biggest discrepancy occurs 
around logi = —2.5 and is at worst a factor of two. 

2.2. Parameterizing the force multiplier 

Although the force multiplier M can be directly used to 
specify the line force, it is convenient to parametrize its 
dependence on t for use in wind calculations. Although 
this approach means that the full complexity of M(t) 
is not captured, it is widely used because of the relative 
ease of manipulating simply-parametrized forms for M(t) 
when deriving wind solutions. 

The basic ansatz under the ICAKI approximation is to 
fit a power-law to the run of M(t) with t, 



M(t) = kr 



(6) 



where a defines the slope and k the amplitude of M at 
t = 1 (i.e., k = M(l)). To capture the flatt ening of M(t) 
for small t, we follow lOwocki et al.l (|1988f ). and modify 
Eq © such that the force multiplier becomes constant 
at low t (as it must in the optically thin limit), 



M(t,0 = *(£)*" 



(1 



c) 1 



- 1 



(7) 



where r max = r] max (£)t. In Eq ([7]) we now explicitly in- 
dicate that M depends on both t and £. Throughout 
this paper, we will choose to describe the influence of £ 
on M via the CAK parameters fc(£) and r] max (^). Allow- 
ing for ^-dependence in these quantities captures the two 
systematic changes in M(t, £) with £: decreasing r? ma x(C) 
with increasing £ allows the turnover in M(t) to shift 
to higher t with increasing £, while a reduction in fc(£) 
at large £ describes the overall decrease in M for larger 
ionization p arameters. 

We follow iStevens k, Kallmanl (119901 ) in choosing that 
a does not vary with £. Although it is certainly pos- 
sible to allow a to vary, this mostly just adds unwar- 
ranted complexity to the parametrization. As is clear 
from Fig. [21 the slope of logM versus logt is not con- 
stant, meaning that a best-fit a is in any case a function 
of the range across which it is fit. Therefore, in all of 
our calculations, we do not fit a but rather fix it to the 
value that is required in order to reproduce the correct 
observed terminal velocity for a single star of the appro- 
priate spectral type in a standard CAK theorjQ. For the 
example star discussed in Section |2~T| this is a = 0.57. 

With a fixed, we derive values of fc(£) and r? m ax(£) by 
fitting our computed M(t) curves to the functional form 
given by Eq ([7]). We restrict this fitting to logt < 0, the 
regime in which M(t) is expected to be dynamically sig- 
nificant. To illustrate the accuracy of this approach, the 

1 The terminal velocity computed in a single star wind calcula- 
tion does also depend on k, but to a much lesser extent than a, as 
one would expect from the functional form of M(t). 



derived fits from our example calculation are over-plotted 
in Fig. [2j As expected, the fits are always very good in 
the optically thin limit and generally agree to within a 
few tens of per cent across the range of interest (i.e. when 
M > 1). However, there are clear imperfections, partic- 
ularly in cases where the slope of M(t) deviates from a 
constant power law (e.g. in our log£ = 2 case). Nev- 
ertheless, the parametrized form provides a convenient 
description and reproduces the force multiplier to within 
a factor of two, which is adequate precision for the pur- 
pose of this investigation. We provide tabulated values 
of fc(£) and /y m ax(£) in the Appendix. 

2.3. Reseating of the force multiplier 

As mentioned in the previous section, in fitting fc(£) 
and ?7max(£) we specified the value of a a priori with the 
aim that the resulting M(t, £) produced a terminal wind 
velocity in agreement with observed values. We now also 
rescale the force multiplier M(i, £) to produce a wind 
mass-loss rate in agreement with observed values, which 
equates to multiplying fc(£) by a correction factoiQ. At 
the cost of some subjective rescaling, this approach has 
the advantage of ensuring that the line force used in the 
colliding-winds model in the following section will pro- 
duce sensible wind parameters while also allowing the 
influence of X-ray irradiation to be explored. We note 
that this modification is of smaller magnitude than the 
current uncertainties in ma ss-loss rates and wind accel- 
eration in massive stars (see lPuls et aLll2008l for a recent 
review) . 

2.4. Results of the line force calculation 

Line force calculations were performed for two massive 
stars: 06V and 04III (see Table Q] for full sets of stellar 
parameters). In Fig. [3] we show the resulting fc(£) and 
?7max(£)- Clearly, for log(£) > the line force is effectively 
suppressed, whereas for log(£) < ionization effects are 
negligible and radiative acceleration is largely unaffected. 

3. THE WIND MODEL 
3.1. Radiatively driven winds 

To compu t e the wind acceleration we follow 
IStevens et alj (|1992f ). Alterations have been made 
to couple the model with a means of estimating the X- 
ray luminosity from the wind-wind collision shocks, and 
then allow for an ionization parameter (£) dependence 
of the line force. The calculation proceeds by solving for 
the wind of one of the stars, with the influence of the 
companion star appearing in the effective gravitational 
potential and as a contribution to the total radiative 
flux in the line force. Subsequently, we change to the 
frame of reference of the companion star and solve for 
its wind in an equivalent manner. In the following we 
describe the solution procedure for each wind. In § 13.21 
we discuss how to infer the shock properties from the 
two wind solutions. 

To simplify the problem we consider steady state so- 
lutions for the flow along the line-of-centres between the 
stars with the forces arising due to orbital motion ig- 

2 The correction factors are 0.74 and 0.34 for the 06V and 04III 
stars, respectively. 
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TABLE 1 

Parameters used for stellar model atmospheres 



Model 


Toff 


Mm 


R* 


log(X/L ) 


log(9) 


Z 


M 




fc(C = o) 


a 




(K) 


(Me) 


(«©) 








(M yr-!) 


( kms" 1 ) 






06V 


38500 


31.7 


10.2 


5.3 


3.92 


1 


2 x 10~ 7 


2530 


0.12 


0.57 


04III 


41500 


48.8 


15.8 


5.8 


3.73 


1 


5 x 10" 6 


2750 


0.18 


0.63 




Fig. 3. — Plots of the line force parameters k (upper) and 
log(f7 max ) (lower) as a function of ionization parameter, £. 



norecQ. We assume that the flow is symmetric about 
the line of centres and that the wind flows purely ra- 
dially from the star. We also assume that the wind is 
isothermal with temperature, T = 0.8 T c ff , where T e g is 
the effective stellar surface temperature. Lastly, we do 
not consider stellar radiation reflected from the oppos- 
ing star's photosphere. Some approximate expressions 
are used in the model to keep the calculations straight- 
forward, whilst achieving an accuracy at the order unity 
level, and in § [7] we discuss possible alternatives. The 
equations for mass and momentum conservation in the 
wind are, 



pv ■ dS = -> Mq = r 2 pv 



(8) 



3 I Parkin et al.l (120111 ) considered the effect of centrifugal accel- 
eration due to orbital motion on the wind acceleration and found 
that it made a correction of a few per cent. Furthermore, for the 
systems considered in this paper, the wind speeds are sufficiently 
large compared to the orbital velocities that there should be no 
significant offsets in th e position of the wind-w ind collision region 
due to orbital motion (Parkin & Pittard 2008). 



F(r, v, dv/dr) 
F(r,v, dv/dr) 



dv 
i — 
dr 



dr 



2a 2 



5rad,(9) 



= 0, 



where r is the distance along the line of centres measured 
from the centre of star 1, v is the wind velocity along the 
line of centres, p is the density, Mq is the mass-loss rate 
per steradian and a is the isothermal speed of sound. 
The gravitational potential due to both stars, 

^ GM^l-T,) GA/» 2 (l-r 2 ) ^ 
r d SG p r 

where M^i and M* 2 are the respective masses of star 1 
and star 2, Ti and r 2 are the respective Eddington ratio 
for each star (Fj = <7 e -£<*i/4'7rGM*ic), d sep is the separa- 
tion of the stars (measured between their centres), and 
G is the gravitational constant. The combined radiative 
line force from both stars, g rac j takes the form, 



#rad = (FlKi 



FzKo 



(11) 



where F\, F 2 are the radiative fluxes, and K%, K 2 are the 
finite disk correction factors (FDCFs) for stars 1 and 2, 
respectively. M(i,£) is the line force multiplier (Eq[7]), 
which in our formulation has a dependence on both op- 
tical depth, t (Eq and the ionization parameter, £ 
(EqEJ). 

The FDCF is a multiplicative factor used to correct 
the point so urce approximation for the finite size of th e 
stellar disk (jCastorllf 97i ICAKl : IPauldrach e Fall [19871 . 
For our adopted geometry and assumptions about the 
flow along the line of centres, we have, 



Ki(r, v, dv/dr) 



(1 + Oi) 



l+Q 



(1 + oinl 



\l+a 



oi(l + a)(l + oO a (l-A&) 



(12) 



where p,\ = cos 9{ with #j being the angle subtended by 
the respective stellar disk viewed from a point in the 



wind, jUjj = 1-R^/r 2 and p'i 2 
o\ = (r / v )(dv / dr) — 1 and a? = 
Following lStevens fc PollockT(fl99 



(d 



~L-RZ2/( d sep-ry, and 



m -r)/v)(dvldr)-l. 
and IPauldrach et al.l 
1986f ) we approximate the FDCFs as purely radial func- 
tions and neglect any velocity or velocity gradient terms. 
In this limit, 



Ki(r, v, dv/dr) — > Ki(r) = 



2 ll+a 



i - [i - /4] 

(l + a)(l-/&)' 



(13) 



In the model considered here, we enforce monotonicity 
in the flow by ensuring dv/dr = m&x(dv/dr,0), which 
ensures that the optical depth parameter, f is a posi- 
tive valued scalar variable. It follows that our models 
do not permit radiative braking ( which requires an in - 
flection in the velocity gradient). iGavlev et al.l ([1997D 
comment that correctly accounting for non-monotonicity 
in the FDCF allows radiative braking. Another way of 



6 



E. R. Parkin & S. A. Sim 



viewing this is that radiative braking involves a bridging 
between two m onotonic flows, which is not facilitated by 
standard C AK theory. 

To proceed, we make a coordinate tra nsform using the 
substitution of variables (jAbbottlll980( ). 



-2GM 1 (l-r 1 ) 



T^^GM^i-rj]- 1 , 

dr 



(14) 
(15) 
(16) 



leading to 
F(u, w, w') 

where 



w' + h(u) - gE{w', £,)B(u)w' a , 

(17) 

(18) 



E(w',£) = k(£)C- a w' 
C = 

B(u) = 
h(u) 



, 9 = r 1 /(i-r 1 ), 
l 



?7maxC 

GM. 1 (i-r 1 )' 
/M* 2 r 2 



l-Q 



77 - <Jh 



1 + 1 

u 



\M ifl Y 1 
M* 2 (l 



X 2 (u)A(u) 

r 2 ) 



and 



A(u) 



A(«) 



(19) 
(20) 

(21) 

(22) 

(23) 



(Note that our definition of C in Eq (|2"0|) differs from 
iStevens fc Pollockl (fl99l 's equation (13)). The FDCFs 



Ki(u) = 



1 - [1 - {u/u^f 



and 



K 2 (u) 



(I + a)(u/u^) 2 ' 

1- [l-(u/u* 2 ) 2 A(u)] 1+a 
{l + a){u/u* 2 ) 2 A{u) ' 



(24) 



(25) 



where = u(d scp ), tt*i = u(i?*i). The strategy for find- 
ing a consistent wind solution is centered around the use 
of the critical point conditions, 



f 1 (w c ,w' c ,C c )=F(u,w,w') =0, 
OF 

f 2 (w c ,w c ,C c ) = — = 0, 



9w 



0. 



(26) 
(27) 

(28) 



Eqs (f26|) - (|28|) are, respectively, the equation of motion, 
the singularity condition, and the regularity condition. 
The subscript "c" denotes the value of the given param- 
eter at the critical p oint. Our set of equatio ns differs 
slightly from those of IStevens fc Pollockl (|1994D . Specif- 
ically, we lack the singular presence of the eigenvalue 



of the problem (namely the mass-loss rate). There- 
fore, we cannot use the equations for fi, f 2 , and /3 
to derive closed form relations for w c and w' c . Instead, 
noting that Eqs (f2"6"| - (|2"5|) compose three equations in 
three unknowns, we solve for w Cl w' c , an d C c using 
a mu lti-dimensional root finder (see, e.g., I Press et al.l 
19861). where the critica l point conditions derived by 



Stevens fc Pollockl (fl99l are used as the initial guess. 



3.2. The post-shock winds 

Once the wind profiles have been calculated we pro- 
ceed to estimate the X-ray luminosity from the individ- 
ual wind collision shocks. The separate values are then 
combined to evaluate the total shocked-wind X-ray lumi- 
nosity. The final step is to use the estimate of the intrin- 
sic X-ray luminosity from the shocked winds fEq I29[) to 
evaluate the ionization parameter, £ (EqEJ). (Note that 
when we swap from the frame of reference of one of the 
stars to its companion's, wc interchange the indices in 
Eqs [21 and Gp. 

We approximate each shock as a thin shell, and take 
, the pre-shock wind velocity and density to be w s h = 
v(r h& \) and p sh = p(r b ai), respectively, where r ba i is the 
ram pressure balance point. The mean post-shock gas 
temperature (i.e. averaged over the bow shock), T ps , 
can be estimated from the Rankine-Hugoniot shock jump 
conditions, kT ~ | X 1.17(w s h/10 8 cm s -1 ) 2 keV, where 
the factor of a half is a correction to account for shock 
obliquity. The 0.01-10 keV X-ray luminosity from the 
respective wind-wind collision shocks is then estimated 
using the simple relation: 



1 • 2 



1 



1 + X 



(29) 



where the total mass-loss rate is approximated as M w 
AttMq. (Note that the X-ray emitting region of the 
shocks is taken to be a point source situated at the ram 
pressure balance point along the line-of-centres) . The 
parameter 5 approximates the thermalization of wind 
kinetic power. In § [4] we consider models of wind-wind 
collision and wind-photosphere collision (as a result of 
the stronger wind overwhelming the weaker wind). In 
the latter circumstance we take S to be the fractional 
solid angle subtended by the disk of the companion star, 



1, 

5 =-(l-c. 



(30) 



where e 2 = 1 — (R« 2 / d scp ) 2 . For a wind- wind collision we 
take S to be the fractional wind kinetic power normal to 
the contact discontinuity ([Zabalza et al.ll2011[ ). 



1 



4 Vl + Coff 



(31) 



where the effective wind momentum ratio of the system, 



off 



M 2 v sh2 

M!V shl 



(32) 



The cooling parameter, x appearing in Eq (|29|) derives 
from the ratio of the characteristic flow time to the cool- 
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Fig. 4. — The number of iterations, iViter required to reach con- 
vergence in the O6V+O6V binary model calculations. 

ing time (|Stevens et al.lll992t ). 

- ( Vsh V ( dsc P \ ( M ^ 
X VlO^ins" 1 / V 1Ql2 cm / ylO- 7 M Q yr- 1 J 

(33) 

If X < 1 the post-shock gas is radiative, whereas if x ^ 
1 the post-shock gas is adiabatic. It is useful to note 
the different scalings of Lx as the importance of cooling 
changes. For adiabatic shocks we have Lx oc (M/v s h) 2 , 
therefore a decrease in w s h leads to an increase in Lx- In 
contrast, when the shocks arc radiative Lx oc Mv^ h , and 
a decrease in v s h reduces Lx- 

In § U we consider calculations with either an attenu- 
ated or unattenuated X-ray flux. For the latter we merely 
have, F x (r) = (L X i +ix2)/47r(r ba i -r) 2 . For the former, 
an attenuated flux is calculated by first scaling a 0.01-10 
keV X-ray spectrum - derived from th e MEKAL plasma 
code (|Kaastralll992l: IMewe et alll995ft - such that its to- 
tal luminosity matches the value from Eq (|2l?)) . The spec- 
tra from both winds are then combined, and the column 
density of gas upstream of the shock, Nn = J rbal pir)dr, 
is used to attenuate the spectrum. (Absorption due to 
the post-shock layers is neglected.) Finally, the resulting 
spectrum is integrated to acquire the attenuated lumi- 
nosity, ixatt from which the ionization parameter can be 
determined. T o this en d we use version c08 .00 of Cloudy 
()Ferlandll2000l see also lFerland et al.lll998h to calculate 
the opacity. 

To test the accuracy of our model, we made a calcula- 
tion for an 06V+06V binary at a separation of 30 R©and 
com pared the estim ated X-ray luminosity t o mode l cwbl 
from|Pjttard| (120091 ) and lPittard fc Parkin! (|2010D (a 3D 
hydrodynamical model with radiatively driven winds). 
We found that our model over-predicted the intrinsic 
0.1-10 keV X-ray luminosity by a factor of roughly two. 
Therefore, for all models examined in this paper we mul- 
tiply Eq (|29|) by a factor of 1 /2. Furthermore, when com- 
puting log(Lx/i , boi) from our models we define Lx as the 
0.5-10 keV X-ray luminosity to be consistent with obser- 
vational studies (see, for example, iNaz c 2009; Naz e et ah! 
[20111: IGagne^Fal1[20rl . 

3.3. Solution strategy 
To summarise, the steps in the calculation are: 



10 15 20 25 30 

r (Rsun) 

Fig. 5. — Wind velocity as a function of radius for the 06V binary 
with a separation of d Bcp = 60 Rq . The different curves show the 
consecutive iterations. Convergence is reached after ~ 5 iterations. 

1. Begin by setting £ = everywhere. 

2. Compute v(r) and p(r) using the wind acceleration 
model described in § 13.11 for the respective stars. 

3. Determine the ram pressure balance point between 
the winds and use this to find Lx for both post- 
shock winds. For a wind-photosphcrc collision, the 
balance point is taken to be at the surface of the 
companion star. 

4. Calculate the X-ray ionization parameter, £(r) (see 
§& 

5. Calculate the change in u s h relative to the last it- 
eration and if convergence is not achievecQ then 
repeat steps ©-©• 

Fig. |4] shows the number of iterations required to reach 
convergence. A larger number of iterations are required 
for smaller separations where the affect of SRSs is great- 
est. 

4. RESULTS 

In this section we examine the influence of the ion- 
izing X-rays from the wind-wind collision shocks on the 
resulting wind acceleration. We have constructed two bi- 
nary systems which wc use to explore the impact of self- 
regulating shocks across a small range of spectral types: 
an 06V+06V binary and an 04III+O6V binary. We also 
consider the collision of the 04III star's wind against the 
photosphere of the 06V star, which arises at binary sep- 
arations smaller than 300 Rq in the 04III+06V case. 
In the following sections we first examine the general 
properties of the self-regulating shocks scenario using the 
06V+06V binary system as our fiducial test case, then 
consider the importance of SRSs as a function of stellar 
separation for the different model binaries. 

4.1. General properties 

We begin by examining the 06V+06V binary at a 
separation of d scp = 60 Rq. Fig. [S] shows the dramatic 

4 For the calculations presented in this paper we required the 
fractional difference in v s ^ (summed over both winds) between con- 
secutive iterations to be < 10 -4 . 
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Fig. 6. — Radial profiles of £ (upper), k(r) (middle), and the 
total column density measured from the shock through the wind, 
Nn (lower). 



Fig. 7. — Comparison of calculations for the O6V+O6V binary 
with and without self regulating shocks. Orbital periods are calcu- 
lated assuming circular orbits. From top to bottom: M, v^, and 
X- 



influence of SRSs on the wind acceleration compared to 
a calculation without this effect included. Three results 
are immediately apparent from this plot: i) the pre-shock 
velocity is considerably reduced, ii) the wind acceleration 
is inhibited, and, iii) the acceleration region is smaller. 

What causes such a significant difference between wind 
calculations with and without SRSs? Fig. [6] shows the 
variation of £ with radius from the star. Close to the 
shocks (which reside at r = 30 Rq in this example), 
log(£) ^> which is sufficient to strongly suppress the 
line force ( § EJ) - In fact, throughout most of the wind 
the value of £ is large enough that the wind acceleration, 
ffrad fx will be inhibited somewhat (see Fig- |3j) , and 
this becomes clear when one examines the run of fc(£) 



against radius. We note, however, that the sharp rise 
in £ close to the shocks is an unphysical consequence of 
concentrating all of the X-ray emission from post-shock 
winds at the stagnation point. SRSs are important when 
£ reaches relatively high values in the wind acceleration 
region (i.e. well away from the shocks) and, therefore, 
the spike seen in the top panel of Fig. [5] does not affect 
the main conclusions of this work. 

For comparison, we also show in Fig. [6] (top panel) a 
calculation performed with an unattenuated X-ray flux, 
which differs only very slightly from the calculation with 
an attenuated X-ray flux. Although the total accrued 
column density, -/Vntot steadily increases when tracking 



Self regulated shocks in massive binaries 



9 



, (days) 



10 



> 

CD 

i 1 



10 



100 



1 1 1 1 1 1 1 


1 1 1 1 . 








/ SRS 




/ no SRS 




1 Observations + 




-7.5 



1000 



d se „ (Rsun) 



Fig. 8. — Plots of kT (upper) and log(Lx/£bol) (lower) against 
binary separation and orbital period (assuming circular orbits) for 
the O6V+O6V binary system. Details of the observations are given 
in §|6] 



back from the shocks towards the star (lower panel of 
Fig. it never reaches a sufficiently high value to im- 
pact the X-ray flux from the wind-wind collision shocks. 
Therefore, the decrease in £ as one moves away from the 
shocks results from an increase in the wind density and 
geometrical dilution of the X-rays. This result holds true 
for all of the models considered in this work. However, 
the influence of attenuation may become more important 
for higher mass- loss rates and/or when the winds contain 
optically thick clumps. 

The wind mass-loss rate is largely unaffected by SRSs. 
This is because the mass-loss rate is set very close to the 
star and, as is evident from Fig. El the ionization parame- 
ter is low enough in this region to have little affect on the 
line force. It is interesting to note that the influence of 
SRSs inhibits the wind acceleration which causes the in- 
ner wind density to increase, thus reducing the influence 
of X-ray ionization (£ oc p^ 1 ). 

4.2. Variation with binary separation 

To better understand the region of parameter space in 
which SRSs will influence the dynamics of the flow and 
the observable properties of a binary system, we have 
performed further model calculations for the 06V+06V 
binary at a range of binary separations, the results of 
which are shown in Figs. [7] and [5J As anticipated from 
§ 14. H and illustrated by the calculations with and with- 
out SRS, the mass-loss rate does not vary greatly due 



to SRSs. The decrease in mass-loss rate with decreas- 
ing binary separation occurs due to the inhibition of 
wind acceleration by the opposing star's radiation field 
(IStevens fc Pollocklll99lFI . 

Also evident from Fig. [7] is that SRSs reduce the pre- 
shock wind velocity for a large range of binary separa- 
tions - for the 06V+06V binary the relative difference 
in w s h between calculations with and without SRSs is 
45% at d sep = 40 R Q and steadily decreases to 6% at 
rfscp = 1000 Rq. SRSs also cause the downturn in u sn to 
occur at larger separations than without SRSs. 

A secondary effect of a lower v s h is that the impor- 
tance of radiative cooling increases (lower \)- As Fig. [7] 
illustrates, without SRSs we would not expect radiative 
shocks (x 1) until d sep < 40 R (P orb < 4 days). In 
contrast, with SRSs this range increases out to d scp < 
70 R Q (Porb iS 10 days). Therefore, we expect that 
systems will have radiative shocks for larger separations 
than previously anticipated. 

The inclusion of SRSs causes a reduction in plasma 
temperature, kT, at all separations but only significantly 
affects the X-ray luminosity, Lx, for a limited range 
of separations (Fig. [8|). At large separations (d scp > 
200 R©) we do not predict a considerable difference in 
Lx due to SRSs. Within our model, this stems from the 
scaling Lx oc (M /w s h) 2 for adiabatic shocks (see Eqs (|29|) 
and (|33| ). and the differences in M and v s h between mod- 
els with and without SRSs (Fig. [7J . These competing ef- 
fects effectively cancel to produce almost identical X-ray 
luminosities. For example, for separations greater than 
200 R© there are uniform offsets of roughly 7% for M and 
v s h between calculations with or without SRSs. 

There are, however, a range of separations where SRSs 
are predicted to make the system intrinsically brighter. 
For the 06V+06V binary this range is 70 < d scp < 
200 R©, and arises because v s h decreases more rapidly 
with decreasing binary separation with SRSs than with- 
out (and because Lx oc (M/v s h) 2 at the relevant values 
of x) ■ With SRSs, and at d sep < 70 R© , \ < 1 therefore 
Lx oc Mv 2 h and the decrease in w s h caused by SRSs re- 
duces the intrinsic brightness of the wind-wind collision 
shocks. The abrupt flattening of Lx as the separation is 
reduced below d ~ 70 R in the model with SRSs (Fig. [8]) 
is a consequence of the transition from \ > 1 to x < 1 ; 
which alters the dependence of Lx on v s h. 

4.3. 04III + 06V binary 
We now consider an 04III+06V binary with unequal 
wind momenta. Based on values for M and Voa from 
Table [T] (which are calculated using isolated single star 
wind models) the wind-wind momentum ratio, £ = 0.04 
(in favour of the 04III star). Assuming terminal velocity 
winds (i.e. neglecting radiative inhibition and SRSs), the 
distance of the wind-wind momentum balance point from 
the star with the weaker wind is, 



r 2 = d s . 



i + C 



1/2 



scp 



£1/2 



(34) 



Stevens & Pollock (1994) note that for a ratio of Eddington 
factors, ri/T2 J; 4 the mass-loss rate will decrease rather than 
increase compared to the single-star case. Furthermore, from their 
equation (24), it is clear that as d decreases, A c increases, causing 
a reduction in Mq 
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Setting T2 = i?*2 = 10.2 Rq, one estimates the wind- 
wind collision to remain away from the surface of the 
06V star for separations greater than 61 Rq. If we im- 
prove on this estimate using our wind model without the 
inclusion of SRSs we find a stable wind-wind collision 
down to separations of 80 Rq (dashed vertical lines in 
Figs. [9] and [TO]) . However, when SRSs are included a 
stable wind-wind collision is not predicted to occur for 
separations smaller than 300 Rq (dotted vertical lines in 
Figs. 151 and ITU)) . The reason for this drastic increase is 
that SRSs tend to make the weaker wind even weaker as 
it is closer to the source of the X-rays at the wind-wind 
collision. The tendency for SRSs to reduce the strength 



of the weaker wind, therefore, becomes more pronounced 
as the separation of the stars is reduced. Consequently, 
even for comparatively large separations, wind- launching 
fails. This general result states that SRSs will cause a 
wind-photosphere collision in unequal winds systems up 
to larger separations than otherwise expected. 

For sufficiently large separations (d scp > 300 R Q ) a 
wind-wind collision is predicted to occur, and in this 
regime SRSs introduce a minor reduction in M's of 
roughly 4% for both stars (Fig. [5J. SRSs also cause 
a reduction in u sn for both winds, most notably at 
smaller separations, where a sharp downturn arises in 
w s h for the 06V wind, signifying the sudden failing of the 
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Fig. 10. — Plots of kT (upper) and log(Lx/£bol) (lower) against 
binary separation and orbital period (assuming circular orbits) for 
the 04III+06V binary system. The vertical lines indicate the lim- 
iting separation below which a model solution could not be attained 
for cases with (dotted line) and without SRSs (dashed line). De- 
tails of the observations are given in § [6] 



wind-wind collision as the 06V's wind becomes increas- 
ingly weakened by the X-rays from the shocks. Some- 
what surprisingly, although we anticipate that a wind- 
photosphere collision will ensue for relatively large sep- 
arations, the post-shock gas is expected to be adiabatic 
(x ^> 1). This differs from previous models in which a 
wind-photosphere collision is typically accompa nied by 
highly radiative shocks from the weaker wind (|Pittardl 
[1991 IParkin fc Gossedl20Tl . 

Similar to the 06V+06V binary, Lx is largely un- 
affected by SRSs for the 04III+06V binary (Fig. [TO]). 
A noticeable reduction in kT values is, however, intro- 
duced particularly for smaller separations. In this case 
we find that SRSs introduce an offset in kT values for 
separations larger than 7OOR0(P O rb > 300 days), and 
that for closer separations SRSs cause a sharp downturn 
in kT, reflecting the behaviour of v s h for the 06V's wind 
- Fig. [HI Lx remains similar between the models, with 
SRSs causing an upturn in Lx for d scp ~ 300 — 600 R© . 

We note that although the 04III star has the stronger 
wind, the X-ray emission is greater from the shocked 
OGV's wind because a larger fraction of its wind is 
shocked at an angle close to the shock normal (thus con- 
verting a larg er fraction of its kinetic energy into ther- 
mal energy - iPittard fc Stevens] 120021 ). It follows that 
the observed kT will also be predominantly weighted by 
the weaker wind. For the specific parameters used in 
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Fig. 11. — Comparison of calculations for the 04III+06V binary 
with a wind-photosphere collision and with and without self reg- 
ulating shocks. Orbital periods are calculated assuming circular 
orbits. From top to bottom: M, v s ^, and \. (Results are only 
shown for the 04IIFs wind because the 06V's wind is assumed to 
be suppressed.) 

our model, and at separations less than 1000 R Q (P rb < 
400 days) the ratio of X-ray luminosity from the winds 
is 2:1 in favour of the 06V. 

4.4. O4IH + 06V binary with a wind photosphere 
collision 

As mentioned in the precedin g section, the inclusion o f 
SRSs (and radiative inhibition IStevens fc Pollock1ll994| ) 
considerably increases the range of binary separations 
where a wind-photosphere collision will occur in a mas- 
sive star binary system. For our 04III+06V model we 
found that a ram pressure balance between the winds 
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Fig. 12. — Plots of kT (upper) and log(Lx/£bol) (lower) against 
binary separation and orbital period (assuming circular orbits) for 
the 04III+06V binary with a wind-photosphere collision. Details 
of the observations are given in § \E\ 

could not be achieved for separations less than 300 Rq. 
In this section we consider the wind-photosphere collision 
occurring at d sop < 300 R© . For this purpose we use the 
model described in §[3] with the difference that the 06V's 
wind is not included and the shock is assumed to occur 
at the surface of the companion star, r = d scp — R*2- 
The fractional wind kinetic power that is thermalized, 
5 is approximated by the solid angle subtended by the 
06V star as viewed by the 04III star (see § I3.2[) . 

As is clear from the plots of M, u s h, and \ m Fig. [Til 
and kT and log(Lx/£boi) in Fig.[12j SRSs have very little 
affect on the wind-photosphere collision. This is because 
log(£) < in the inner wind acceleration region which 
allows the wind to accelerate to a similar velocity to the 
case with no SRSs. log(£) > is only reached in regions 
beyond the acceleration zone, meaning that the driving 
is unaffected. To illustrate this, in Fig. [13] we show ra- 
dial profiles of wind velocity, log(£), and fc(£) computed 
for binary separations of 50 and 200 R©. Clearly, the 
line force is only suppressed by SRSs in regions where 
log(£) > 0, which reflects the almost step function like 
behaviour of fc(f) at log(£) ~ - 1 (Fig.0). 

5. AN APPROXIMATE INDICATOR FOR SELF 
REGULATING SHOCKS 

It would be useful to have a simple means of estimat- 
ing the separation at which we expect SRSs to play an 
important role in the wind-wind collision. This could be 
used, for example, to estimate whether SRSs should be 
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considered when modelling a specific system. Surveying 
the results of our model calculations, we find that, to 
within an accuracy of a factor of two, we have the fol- 
lowing approximate relations for our equal winds binary 
system: 



V2gm»7r7, 



J- T X "scp 

M 



p(2R*) 



16irv csc Rl ' 



(35) 
(36) 

(37) 
(38) 
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where p(2R*) is the wind density at a radius of 2 i?„ 
and do = 0.5 R© is a constant used to fit the variation 
of S/(l + x) with d. With the further simplification 
of neglecting any attenuation of X-rays as they travel 
back through the wind towards the star (which has been 
shown in § 14.11 to have minor influence, at least for the 
system parameters considered), inserting Eqs (|35[) 
into Eq (|5|) gives, 



£(2i*.) = 647rVm H do^(^-l 



(39) 



From the results of § [4] a requirement for SRSs to affect 
the wind-wind collision is that log(£(2ii*)) > 0. Set- 
ting £(2i£«) = 1 and re-arranging Eq (|3"9")l leads to a 
cubic equation for d sep which, for parameters pertaining 
to our equal winds 06V+06V binary (§§ Qand HH , 
has one real root of d sep ~ 91 R . Inspecting Figs.[7]and 
M one sees that this value is consistent with the onset of 
a marked difference due to SRSs. 

6. COMPARISON AGAINST OBSERVATIONS 

To facilitate a comparison of our model results against 
observed O+O binaries we ha ye extracted a samp le of 
systems from the studies by IGagne et al.l (|201lD and 
iGagne et al.l (J2012) . Only systems with orbital periods 
within the range of our models have been considered. 
We then separated the remaining systems into those with 
roughly equal winds or unequal winds systems, where we 
classify the former as systems in which the stars differ 
by less than a spectral type and/or subclass, and the 
latter as systems which differ by more than this incre- 
ment. The roughly equal winds systems (three in total) 
are then compared to our 06V+06V binary and the un- 
equal winds systems (six in total - although kT values are 
only available for four systems) against the 04III+06V 
binary. 

We remind the reader that the model used in the cur- 
rent investigation includes simplifications to the physics, 
which have been chosen so as to allow a tractable ini- 
tial exploration of the SRS effect and its potential im- 
portance. Nevertheless, in this section we compare our 
model results against observations to provide a sense 
of how SRSs might be relevant in explaining general 
trends. However, we caution that detailed comparison 
must await more thorough modelling. 

6.1. Equal winds systems 

In the equal winds case we find reasonably good agree- 
ment between our 06V+06V model and the observed 
kT and log(Lx/^boi) (Fig. [8]). Including SRSs improves 
the match to the observed kT values. The 06V+06V re- 
sults (Fig. [8} show that for binary separations less than 
~ 200 R© (P rb < 30 days) we expect roughly equal 
winds systems to be brighter than the expected luminos- 
ity from embedded wind shocks in th e resp ective stars 
(log( LxALhoi) ^ -7 , e.g.. iSana efld1l2006b iNaze et al.l 
120 111 ). Roughly equal winds systems with separations 
larger than ~ 200 R© will not, therefore, be identifiable 
as CWB systems from their X-ray luminosity but instead 
they may be identifiable by kT > 0.6 keV (i.e. hotter 
plasma te mperature than anticipat ed from a single mas- 
sive star - lOwocki fc CohenlH999l ). Indeed, it may be 



the case that only early-type O+O binaries with inter- 
mediate orbits and strong winds will be pr olific X-ray 
emitt ers (e.g., Cyg OB#9 - P or b = 858 days - INaze et al.l 
l2012f ). with the majority of later- type massive binaries 
only being identifiable as CWBs (in X-rays) via plasma 
temperatures above 0.6 keV. 

There is significant scatter in the observed kT values 
and log(Lx/^boi) for orbital periods less than six days. 
We do not attempt to compare our model against these 
systems because, when the separation of the stars be- 
comes comparable to their stellar radii, one expects ad- 
ditional effects that we have not considered to become 
important. For example, tidal deformation, gravity dark- 
ening, p hotospheric reflection, and the possibility of mass 
transfer dGavlev et alj|1999t iDessart et aT]|2003b lOwockil 
12001 [Dermine et al.ll2009[ K 

6.2. Unequal winds systems 

At the separations of the observed unequal winds bina- 
ries, a wind-wind collision is predicted from models with- 
out SRSs while a wind-photosphere collision is expected 
based on our SRS calculations. Comparing Figs. [TU] and 
[T2l one sees that both cases do arguably similarly well 
at matching the observations - although the wind-wind 
collision with no SRSs does appear to over-predict the 
observed log(Lx/£boi)- Interpreting this comparison is 
complicated, and it may simply be indicating that reality 
lies between these two different cases, i.e. a wind-wind 
collision prevailing to smaller separations but with some 
wind suppression due to SRSs. 

Examining the wind-photosphere collision in more de- 
tail, the models systematically over-predict kT values by 
roughly a factor of two, irrespective of whether SRSs 
are considered or not (Fig. [T2"|) . However, the agree- 
ment between the model and observations is reasonably 
good for log(Lx/£-boi), with the exception of the two sys- 
tems with orbital per i ods of roughly 6 days: HD93205 
Townslev etlll I2I3T1I: INaze et al.l l20Tl and HD101190 



■Chlebowski et al.1 119891: ISana et alJl20lTt IGagne et al.l 
l2012f ). We remind the reader that orbital periods have 
been converted to binary separations under the ba- 
sic assumption of circular orbits, which is accurate for 
the majority of the systems in the sample. Consider- 
ing the two outliers with orbital periods of roughly 6 
days, the former, HD93205, has an orbital ec centricity of 
0.37 dMorrell et al.ll2001t iRauw et al.l [20091). Using the 
ephemeris from iMorrell et al. I (|2001D and the date of the 
Chandra observation of HD93205, we estimate an orbital 
phase of ~ 0.2. As this is relatively close to periastron, 
we cannot appeal to the larger separation that will oc- 
cur at apastron to improve the match against our model 
results. Adopting the recently derived orbi tal solution 
for H D101190 with an eccentricity of ~ 0.3 (jSana et al.l 
1 2 llf ) does not help the agreement between our model 
and its log^x/^boi) datapoint either. A more detailed 
hydrodynamical model of a wind-photosphere collision is 
warranted to investigate the systematic discrepancy in 
kT values and log(Lx/£boi)- 

7. DISCUSSION 

lOwocki fc Gavlevj (fl99ll and lGavlev et~aTl (fl997h have 
argued that for binary systems where a ram pressure bal- 
ance is not expected to occur, the radiation field of the 
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star with the weaker wind may decelerate the incom- 
ing wind of its stronger companion. We note that the 
high values of £ in the vicinity of the shocks (see Figs. [6] 
and [T3j) raises questions about the ability of radiative 
braking to produce a time-steady wind interaction re- 
gion. For instance, perhaps an incoming flow is initially 
subject to radiative braking, but any shock which subse- 
quently forms (and the associated X-ray flux) will sup- 
press the braking force, leading to the dominant wind 
continuing on its path towards the weaker star's pho- 
tosphere. Then, with the X-ray emitting shocks extin- 
guished - or sufficiently weaker/far enough away from 
the point where radiative braking was originally effec- 
tive - the cycle can repeat. However, more detailed hy- 
drodynamical model is required to properly assess these 
points as it may be the case that close to the photo- 
sphere of a star the gas density will be sufficiently high 
that the ionization parameter will be small (either due 
to an intrinsically dense photosphere, a build-up of gas 
behind the shock, or wind strengths weakened by inhi- 
bition/braking) , in which case radiative braking may 
prevail. Therefore, a complicated, and most likely time- 
dependent, competition between radiative braking and 
SRSs may arise. 

In our model calculations we have examined the in- 
trinsic X-ray luminosity. An important related question 
is how the observed Lx would be affected by shock self- 
regulation? For instance, SRSs reduce the post-shock gas 
temperature and, consequently, the energy of emitted X- 
rays. As the susceptibility of X-rays to absorption in- 
creases at lower energies, the X-ray flux that reaches the 
observer may be comparatively much fainter for systems 
where SRSs are effective. Therefore, although our model 
including SRSs overestimates the observed Lx for binary 
systems with orbital periods of a few days to about the 
same level as model cw bl from lPit~a rd (200(|) (see also 
iPittard fc~ Parkin 20l0j), further work is needed to eval- 
uate how the observed (attenuated) Lx is impacted by 
SRSs. This will be an important point to pursue in fu- 
ture work. 

The model adopted for this investigation features a 
number of approximate relations, adopted to keep the 
calculations relatively simple whilst achieving an order- 
of-unity accurate prediction of the influence of SRSs on 
a wind-wind collision. While these approximations have 
been chosen in order to give a simple exposition of the 
SRS mechanism, it is important to note that alternative 
approximations could have been made, whose respective 
merits should be borne in mind for future investigations. 
Firstly, the approximation used to estimate the half- 
opening angle of the bow shock, 6» half = (7rC c ff )/ (1 + Ceff), 
which features in Eq (|31|) for S does not consider the in- 
fluence of the post-shock wind momentum on the global 
shock geometry. Calculations of shock half-opening an- 
gles which include this addition al momentum flux (e.g. 
iCanto et all 119961: lGavlevl[2009l) find that it widens the 
bow shock, leading to a slightly different, and more accu- 
rate, scaling of #haif with £ e ff- Similarly, a more accurate 
expression for S, in the case of a wind-photosphere col- 
lision (Eql30j) could likely be derived using a glo bal mo- 
mentu m flux approach similar to that adopted bv lGavlevI 
(12001 . Secondly, in using Eq (f5T|) to calculate the thcr- 
malization efficiency it is implicitly assuming that all 
wind kinetic energy normal to the shock is thermalized 



and the effect of shock obliquity is not included (although 
a constant obliquity correction is included when calculat- 
ing the mean plasma temperature). We anticipate that a 
more accurate treatment of shock obliquity would intro- 
duce an order unity correction to the results and could 
improve the agreement with observations. Thirdly, we 
do not include the radiative-driving force arising from 
the stellar radi ation field refl e cted b y the opposing star's 
photosphere. iGavlev et al.l (|1999[ ) examined a similar 
scenario in planar geometry and found tha t the radiative 
inhibition effect ([Stevens fc Pollockl Il994l ) was weaker, 
and the mass-loss rate higher, due to the extra acceler- 
ation force from reflected radiation. To include the re- 
flection effect in a geometry such as illustrated in Fig. Q] 
is not trivial, but would be a worthwhile avenue for fu- 
ture work. Reflection could be particularly important for 
the wind-photosphere collision model as it could enhance 
radiative braking and/or impinge on the wind-bearing 
star's wind acceleration. 

8. CONCLUSIONS 

We have presented steady-state wind models for mas- 
sive star binary systems in which the X-ray emission from 
the wind- wind collision shocks modifies the driving of the 
wind, which we term self regulating shocks (SRSs). To 
this end we include a parameterized radiative line force 
with X-ray ionization dependence (derived from line force 
calculations) in our wind model. Our primary result is 
that X-ray radiation from the shocks is found to inhibit 
the wind acceleration and can lead to lower pre-shock 
velocities, which in turn causes the post-shock plasma 
temperature to decrease. In general, SRSs will alter the 
pre-shock velocity if the ionization parameter log(£) > 
within a radius of 2 We believe the qualitative re- 
sults from this investigation to be robust, but note that 
quantitative estimates made from the model may change 
as more complete physics prescriptions are incorporated 
into new models. Caution should be exercised when ex- 
trapolating the results from the sample calculations in 
this paper to specific systems, as to acquire accurate re- 
sults will require a dedicated analysis. 

Despite the presence of an anticipated feedback loop 
between the shocks and the wind driving, the resulting 
intrinsic X-ray luminosity of the shocks is not strongly 
altered by the inclusion of SRSs. However, although not 
examined in this work, lower plasma temperatures may 
render the X-ray emission more susceptible to absorp- 
tion, which could have an impact on the observed atten- 
uated emission. 

We have presented model results for O6V+O6V and 
04III4-06V binary systems computed for a wide range of 
binary separations. For the O6V+O6V binary the main 
difference introduced by SRSs is the reduction in pre- 
shock velocities described above. For the 04III4-06V 
binary, SRSs greatly increase the separation at which a 
wind-photosphere collision (which occurs when there is 
no ram pressure balance between the winds) from 80 to 
300 R©. Furthermore, close to the shocks, where X-ray 
ionization is greatest, the line force can be completely 
suppressed, and we conjecture that this may render ra- 
diative braking ineffective, or highly time-dependent. 

A comparison of our model results to observations re- 
veals that the inclusion of self-regulated shocks improves 
the agreement for plasma temperatures in roughly equal 
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winds systems. However, irrespective of the inclusion 
of self-regulated shocks we find a systematic offset in 
plasma temperatures for unequal winds systems (which 
we model as a wind-photosphere collision, as expected 
for the range of binary separations probed by observa- 
tions). The models show reasonable agreement with ob- 
servations for log(Lx/iboi)- Unequal winds O+O star 
systems with a wind-wind collision are not expected to 
be brighter than their respective stars in X-rays. Such 
systems are predicted only to have a wind- wind collision 
above some cutoff binary separation because at smaller 
separations SRSs prevent a stable wind-wind ram pres- 
sure balance. For our sample 04III+06V system, this 
cutoff is at a separation of 300 R© (P rb > 70 days). 
However, shorter period systems (separations smaller 
than 130 R Q , P or b < 20 days for our 04III+06V model) 
with a wind-photosphere collision should be noticeably 
bright in X-rays (i.e. log(L x /A>ol) > -7). 
This work is a first attempt at modelling the influ- 



ence of X-ray ionization on wind driving in massive star 
binary systems. In closing we suggest a few possible av- 
enues for future work. Developing more realistic models 
requires multi-dimensionality, with 2D models being the 
logical next step. Furthermore, time-dependent calcu- 
lations would be enlightening as one can envisage that 
oscillatory behaviour may result from perturbations in 
the pre- /post-shock flow, and it will be interesting to ex- 
amine whether SRSs can exp lain flaring in massive sta r 
binary X-ray lightcurves (e.g. iMoffat fc Corcoranl[2009| ) . 
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TABLE 2 



FIT PARAMETERS FOR fc(£) AND ?7max(£) 


log(0 




06V 




04III 




k 


log(»? m ax) 


k 


log(r; max ) 


-5.0 


0.348 


7.035 


0.244 


6.705 


-4.0 


0.349 


7.034 


0.245 


6.703 


-3.0 


0.351 


7.029 


0.247 


6.699 


-2.5 


0.353 


7.023 


0.248 


6.693 


-2.0 


0.356 


7.015 


0.250 


6.685 


-1.5 


0.358 


7.003 


0.252 


6.672 


-1.0 


0.351 


6.992 


0.248 


6.662 


-0.5 


0.321 


6.992 


0.223 


6.675 


0.0 


0.263 


6.999 


0.176 


6.722 


0.5 


0.196 


6.966 


0.128 


6.721 


1.0 


0.130 


6.526 


0.090 


6.288 


1.5 


0.103 


4.891 


0.087 


4.268 


2.0 


0.035 


4.164 


0.029 


3.392 


2.5 


0.021 


2.660 


0.018 


2.305 


3.0 


0.023 


2.068 


0.020 


1.911 


3.5 


0.020 


1.655 


0.019 


1.563 


4.0 


0.014 


1.233 


0.015 


1.220 



Note. — The lCAKl parameter a is fixed for each star, respectively (see Table [T] and §[2J. 
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Stevens, I. R., & Pollock, A. M. T. 1994, MNRAS, 269, 226 Zabalza, V., Bosch-Ramon, V., & Paredes, J. M. 2011, ApJ, 743, 
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APPENDIX 

A. TABLES OF FITS TO LINE FORCE PARAMETERS 

In § [5] we described the dependence of the line force on the ionization parameter, £. To allow a straightforward 
application to the wind model in § [3] we described this dependence in terms of the parameters fc(£) and ?7max(£)- 
Tabulated values of these parameters are provided in Table [2] Note that these values have not been rescaled (as 
described in § I2.3|) . For log(£) > 4 the line force is strongly suppressed (M ~ 0). 



